Semiclassical eigenstates of four-sublattice antiferromagnets 
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Applying Bohr-Sommerfeld quantization and the topological phase of spin path integrals, one can 
determine the multiplicities, lattice symmetries, and eigenvalue clustering pattern of the low- lying 
singlet eigenstates of the triangular and fee antiferromagnets with 4-sublattice classical ground 
states. In the triangular case, the clustering pattern agrees with numerical results of Lecheminant 
et al (Phys. Rev. B 52, 6647 (1995)). 
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How can one identify the long-range order of an ex- 
tended quantum system from exact diagonalizations that 
are severely constrained by finite-size effects? The com- 
monest analysis is to extrapolate the ground-state corre- 
lations and energy gaps infinite size. An alternative ap- 
proach uses the finite-size splittings and symmetries of an 
entire family of low-lying states to characterize possible 
symmetry breaking. ^],^[ In the latter spirit, we present 
a semiclassical calculation of the low-lying singlet eigen- 
states and eigenenergies of four-sublattice Heisenberg an- 
tiferromagnets, assuming long-range order in the N — > oo 
limit (where N is the number of spins). The Hamiltonian 



#({*}) = 5 E J ( r < 



(1) 



where spin Sj has quantum spin length s, and ry is the 
vector connecting sites i and j. By "four-sublattice" an- 
tiferromagnet, we mean one for which a classical ground 
state is any state with spins parallel within each of the 
four (equivalent) sublattices, and the vector sum of their 
sublattice magnetizations zero. Important examples are 
the triangular antiferromagnet (|J§], with J\ > and 
second- neighbor coupling J2 G (J1/8, Ji), and the Type 
I fee antiferromagnet [Qj5| with J\ > 0, J2 < 0. After al- 
lowing for rotations, there is still a two-parameter family 
of classical ground states. 

These special degeneracies will get split in a fashion 
characteristic of four-sublattice order, producing a char- 
acteristic pattern of low-lying eigenstates. To derive this 
pattern, we will map the system (approximately) into 
a 4-spin one, and in turn to a sort of one-spin system, 
very much like the well-known cubic-symmetry molecu- 
lar rotor HQ. Using Bohr-Sommerfeld (BS) quantiza- 
tion of classical orbits to define energies, which are split 
by tunneling between the classical orbits, we predict be- 
low a characteristic pattern of level clustering in accord 
with the numerical data Q for the triangular case with 
N = 28, the largest size studied to date. 

Let C a be the set of N = N/4 sites in sublattice a 
(a = 1,2,3,4). The sublattice magnetizations are four 
big spins S a = ^2 ieC s,-. Inspired by the classical ground 
states, we will define "four-spin " states as those in which 
(i) each S a has maximum total spin length |S Q | = S = 



Ns; (ii) the total (quantum) spin S tot = J2 a is zer0 
(singlet). This is the low-energy singlet subspace, as first 
recognized by Lecheminant et al It has 25+1 

independent quantum states. 

Our first goal is a reduced Hamiltonian acting only on 
four-spin states. We first obtain an infinite-range version 
of (^) - i.e., every spin interacts equally with every spin 
of a different sublattice - by replacing each Sj — + S a /N 
in 7i, where i € C a : 
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(2) 



where J = {3N)- 1 £f ff J(r), J = J - 
N-iJ2l ame J(r). (Here "same" and "diff" mean the 
sums are restricted to interactions connecting the same 
or different sublattices.) In the triangular case J = Jq — 
2( Ji + J 2 )/A>. 

The "four-spin " states (having S tot = and maxi- 
mum |S Q |) are manifestly the only ground states of 7ii„/, 
having energy E in f — ~2JS(S + 1). Their degeneracy 
is broken by an effective Hamiltonian Ti. se i({S a }) which 
"selects" particular ground states; this will be crudely 
approximated by a phenomenological biquadratic form 



a<f3 



S/3) 2 + Cf 



biq 



(3) 



with K > 0. 

Eq. (H) can derived from 5Tt = H — Ttmf by either of 
two perturbative approaches. First, in the large-s limit, 
and after the usual Holstein-Primakoff expansion around 
a four-spin coherent state, SH. could be approximated by 
a harmonic Hamiltonian which is diagonalized by spin- 
wave states. The spin-wave zero-point energy, summed 
over all modes §,|l0|, can then be expressed as a func- 
tion of the four spin directions. By general ("order due 
to disorder") arguments, collinear (parallel or antipar- 
allel) spin configurations have the lowest zero-point en- 
ergy H|lO| , and E q. (|3| ) is the simplest analytic form with 
this property Q[ljF 

Alternatively, Eq. (^) can be directly derived from 
second-order perturbation theory in STi.. This approx- 
imation |l3| (related to that of |ll]]) gives precisely the 
form (j^), with 



1 



K = 



J? + Jj-2N- 1 {J 1 + J 2 ) 2 
NJS(2S-l) 2 



(4) 



for the triangular case, and Cu q = — 25 2 (55 2 — 65 — 
2)K. We caution that the (presumably large) quantum 
fluctuations within each sublattice ought to renormalize 
(Q) by a substantial factor. 

From here on we work with the effective four- 
spin Hamiltonian Ei n f +7ibiq- The quantum eigenstates 
of Tdbiq will be approached semiclassically; that is valid 
in a moderately large system even for s = 1/2, since 
5 = Ns is the large parameter in Tibiq. Our classical 
coordinates are the unit vectors m Q that parametrize 
the four-spin coherent states, constrained (just like four- 
sublattice classical ground states) only by J2 a ™« = 0- 
It will be useful to define coherent states |{m Q }) for the 
four-spin problem, meaning each S a has maximum pro- 
jection in the direction m Q . There are competing pre- 
scriptions to derive the effective classical Hamiltonian 
Ubiq from Ti-biq, which disagree after the leading power 
in 5; we adopt @,|| 



({m a }\H Mq \{m a }} = U biq ({m a }) + C L 



where 



U biq ({rh a }) = -K\ 



a</3 



(5) 



(6) 



with K v = 5 2 (5 - l/2) 2 K and C v = 5 2 (65 + b/2)K + 

Cbiq- 

Now, the Green's function of such a system can be 
expressed as a path integral, |l4|,[l5| in which each path 
is weighted as usual by 



exp 



J dtU blq ({m a }) - $ 



where the topological phase of the path is given by 



$ = 



4 
a=l 



Sfi({m Q (f)}), 



(7) 



(8) 



where f2({m(t)}) means the spherical area the trajectory 
rii(f) has swept out on the unit sphere about its "north 
pole" p™. The classical dynamics is 



dm a /dt = yrh a x SU biq ({m a })/dm a 



(9) 



The dynamics in the 5-dimensional subspace ^ m a = 
is separable by a change of variables from {rh Q } to a 
unit vector n and a proper rotation matrix 1Z defined as 
follows. Let n M = i(rh M +rii4), fi = 1, 2, 3 ||; these three 
vectors are orthogonal as follows from ^ m Q = 0. The 
proper rotation matrix 1Z is defined to align these vectors 
(in either sense) along the respective coordinate axes e M , 
so n M = n M e M and n = (m, 112,113) is a unit vector. A 



discrete redundancy remains, that 1Z is well defined only 
up to a 7r rotation of n about any coordinate axis. 

We are finally interested in the singlet projection of 
the coherent states, which is equivalent to averaging over 
all 1Z (with correct phase factors). This singlet basis is 
labeled only by n; it is important for the sequel that 
n's related by the redundancy correspond to identical 
singlet-projected basis states. 

Substituting from (JsJ) shows the Hamiltonian depends 
only on n: 



U biq (n) 



-8K u y(n^) 4 + 2K l 



(10) 



When (j|) is translated into the coordinates n and 7Z, it 
turns out that dlZ/dt = while dh/dt = 7'n x SU b i q /Sh, 
identical to the classical dynamics of one spin with 
Hamiltonian (|i"fj|), i.e., a cubic anisotropy field; these clas- 
sical orbits follow contours of constant energy on the unit 
n sphere, in the sense indicated by the arrows in Fig. [l]. 
Furthermore, each rh a (t) traces out the same shaped tra- 
jectory and makes the same contribution to the sum (^) . 
Hence, the total topological phase $ = 45f2({n(t)}) is 
the same as that of one spin with length 5 = 45 = Ns. 
In (only) this sense, we have mapped a four-spin to a 
one-spin problem. 

Bohr-Sommerfeld orbits and tunnel splittings — We 
can identify three kinds of classical ground state with 
high symmetry, corresponding to the special points in- 
dicated on the n sphere in Fig. |l|: (i) The "collinear" 
(C) states, with all spin directions m a oriented along the 
same direction (two being parallel to it, and the other two 
antiparallel) ; these are the ground states of U b i q . There 
are three C states (since there are three ways to group 
the spins {m Q } into pairs). On the unit n sphere, Cx,y,z 
lie along the ±£, ±y and ±z coordinate axes; these six 
points really correspond to three distinct states in view 
of the discrete redundancy, (ii) "Tetrahedral" (T) states 
in which the spins point towards the corners of a regular 
tetrahedron in spin space, the maximum energy states 
of Ubiq. There are two T states since there is a right- 
handed and a left-handed way of orienting the tetrahe- 
dron. These correspond to points (±1, ±1, ±l)/\/3 on 
the n sphere, corresponding to only two distinct states 
after applying the discrete degeneracy, (iii) "Saddle" or 
"square" (5) states, in which the four spin directions lie 
in the same plane in spin space and differ by 90° rota- 
tions; these are saddle-points of the Uuq function. The 
corresponding values of Uu q (using (^) are Uc — —6Ku, 
U s = -2K V , and U T = -\K V . 

The BS quantization condition W\ selects orbits 
around T with energies Uu q (n) = Uti such that $ = 2nl 
for I — 0,1,2,.... However, the BS condition for the 
energy Uci of a C type orbit is that the complete loop 
has phase $ = Airl; this is due to the redundancy of ii, 
whereby a C type orbit returns to an equivalent state 
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(and completes the true orbit) after looping just halfway 
around the C point. We will label each orbit by its 
BS quantum number I = 0, 1, . . . and by the label of the 
fixed point it encircles, thus "T; 7 " (for 7 = ±, a twofold 
degeneracy) or "C; 7 " (for 7 = X, Y, Z, a threefold degen- 
eracy.) See Fig. |l|. 

Each level cluster is built from BS orbits degenerate in 
energy. Depending upon the topological phase ]l7|], tun- 
neling may occur between these orbits and slightly split 
this degeneracy. Let's write, e.g., tci(Y[c]X) to mean the 
total amplitude to "hop" from orbit Cix to Ciy along 
paths which separate from Cix, pass near saddle point 
S c , and join onto Ciy- (In light of the discrete redun- 
dancy of n, there was just one tunneling path connecting 
a pair of C orbits.) All symmetry-related hoppings have 
the same magnitude tci but their phases depend on the 
"gauge" choice used in defining coherent states. 

The hopping between the three BS states is just like 
that between three atomic orbitals in a ring threaded 
by a flux. The eigenenergies are Uci + 2\tci\ cos((27rj — 
Re®ci)/3), where j is any integer, and the gauge- 
invariant Re$c"i is the cyclic sum of the three phase an- 
gles. We will write "(^1, ^2)" for the pattern of eigenval- 
ues in this cluster, meaning the (lower, higher) eigenvalues 
have degeneracies (^1,^2)- respectively. 

So we just need 



to know t cl e lRc ^' = tci(Z[a]Y)tci(Y[c]X)t C i(X[b}Z) 



the amplitude for closed paths around the loop shown in 
Fig. 1(b). The stationary-phase trajectory, as shown in 
Fig. 1(b), follows orbit Cix to a point along the "equa- 
tor" , then crosses the classically forbidden barrier by fol- 
lowing a trajectory with part complex coordinates and 
joins onto orbit Ciy such that Uuq = Uci along the 
entire path JTBfl . Because the classically forbidden seg- 
ments follow high symmetry lines, the real part of the 
spherical area enclosed (hence of $cz) is computed from 
the real part of the trajectory (shaded in the figure), 
and correspondingly for the imaginary parts. For exam- 
ple, the segment of stationary-phase path crossing S c in 
Fig. [l] has 9{4>) = n/2 — Hm6(<f)); in the spherical area 
il = § d(f>(l — cos #(</>)), the integrand for this segment 
has cos 6* — + isinh(Im#). 

Thus Re$(7i is S times the spherical area enclosed by 
the loop shown in Fig. 0(b): i.e. S , (4tt/8) - 3(4ttZ/4), 
where the second term comes from the arcs truncat- 
ing the triangle corners in Fig. [l] (b). Thus e lRc * C! = 
(— l) s / 2_3i . Recalling that S = Ns, we obtain a cluster 
pattern (2, 1) if Ns/2 — 3/ is even or (1,2) if it is odd. 
On the other hand, Im$c7 is just the WKB exponent ap- 
pearing in the tunnel amplitude, thus tci ~ e~ Im * Ci / 3 . 

The case of T type orbits differs in that there are 
three distinct stationary-phase paths connecting orbits 
Ti + and T;_. So, from the truncated-square loop shown 



[&]... [c] — > [c]...[a] — > [a] ...[&] as well, which give all the rel- 
ative phases between the three different paths connecting 
orbit T; + to T;_. The total hopping is a sum over these 
paths, £y* = J2 P =a b c tTi(— [p]+)j an d the relative phases 
of these terms turn out to be 0, ±Re$rz- So when Ns — 41 
is divisible by 3, then t^f ^ and the cluster pattern is 
(1,1) with splitting 2|i£?/|; otherwise, = and the 
cluster pattern is (2) (i.e., unsplit) pl| . 

The overall clustering pattern thus begins (starting 
with the ground state) with repeats of (2, 1)(1, 2) . . . for 
C orbits - for odd Ns/2, the first cluster is (1, 2) - and 
turns into repeats of . . . 2(1, 1)2 . . . for the T orbits corre- 
sponding to the higher energies among the low-lying sin- 
glets. Amusingly, the sequence of degeneracies {2, 1, 1, 2} 
continues unbroken past the saddle-point energy Us, so 
there is no sharp boundary between the two behaviors. 

Comparison to diagonalizations — The spin-1/2 trian- 
gular system was exactly diagonalized by Ref . (|] . Table | 
shows the numerical data of Lecheminant et al for 
the s = 1/2 triangular antiferromagnet with (J\, J2) = 
(1,0.7), for N = 16 and N = 28. All of their low-lying 
singlet energies are given, compared with our numerical 
predictions, as AE (the difference E — Ei n f). Each row 
is one cluster, labeled in column 1 with the BS orbit from 
which it derives, and in parentheses the degeneracies of 
the levels in the cluster (from lowest to highest). The 
other columns give the mean energy of the cluster, and in 
brackets its tunnel splitting tci or $n (if nonzero). The 
"4-spin" column is (0) with K given by (|j); the "Theory" 
column is (|5|) and ( jl(i| ) with Kjj given after (^J). (In the 
absence of a theory for their prefactor, we estimated the 
"theory" splittings as tci — ► (8Ku)exp(—S Im$ci/3), 
where Im$ c; /3 = 0.55 for I = and 0.10 for I = 1.) 

The eigenvalue clusters (which comprise all the low- 
energy states in these small systems) fall in exactly the 
pattern we predict. For N = 16, there is a near- 
cancellation in formula (0); this explains qualitatively 
why the whole energy scale is so small compared to 
N = 28, and why quantitative agreement between theory 
and experiment is not expected in the N = 16 case. For 
N = 28, the cluster energies agree fairly well, apart from 
a constant offset which we cannot at present explain. 

Summary and discussion - To conclude, we have iden- 
tified two small energy scales among the low-lying singlet 
states in four-sublattice antiferromagnets, also indexing 
all of these eigenstates and explaining the observed pat- 
tern § of their energy splittings. The smallest splittings 
(of 0(exp(— constTVs))) are explained by tunneling be- 
tween different classical wells that result from discrete 
symmetry breakings; exact degeneracies occur when tun- 
nel amplitudes summed along multiple paths cancel, ow- 



in Fig. 0(c), we find t T i(+[b}-) 2 t T i(-[c} 
with Re<I> T; = S , (4tt/6) - 4(2?rZ/3). From the other 
loops related to it by symmetry, we get formulas with 



ing to the topological phase. |17| The next smallest split- 
tings between singlets in the N < 28 systems are between 
clusters, i.e. between successive Bohr-Sommerfcld orbits. 
This spacing is proportional to the zero-point energy per 
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spin (favoring spin collinearity) which scales as a constant 
as N — > oo. 

We can, of course, predict cases which have not yet 
been diagonalized: e.g., for s = 1/2 in iV = 32 
(fee case) or N = 36 (next larger triangular system), 
the cluster patterns are respectively (1, 2)(2, 1)1(2) and 
(2,1)(1,2)(2)(1,1). 

For sufficiently large N, our assumption that spins in 
each sublattice stay rigidly aligned must fail during the 
tunnel event: the tunneling barrier will be smallest for 
an inhomogeneous quantum-nucleates and then grows. 

Finally, it may be noted that our semiclassical treat- 
ment is exactly the same as that of a one-spin Hamil- 
tonian with total angular momentum S and cubic 
anisotropy - exactly the Hamiltonian analyzed semiclas- 
sically HQ to account for the rotational spectroscopy of 
SFg. The eigenstate clusters in Table | are a subset of 
those in the one-spin system, as 3/4 of the latter have 
symmetries forbidden in our system (in view of the dis- 
crete redundancy of states labeled by n). 

We are grateful to P. Houle, K. Nakamura, and 
E. Heller for discussions, and to the authors of Ref. [2] for 
comments and for supplying their numerical data. This 
work was supported by NSF grant DMR-9612304. 



state of any term in S7i, say (s, — S a /N) ■ (sj — Sp/N) 
with i G C a , j £ jCp (with a ^ /?). It leaves S to t = but 
reduces |S a | and \Sp\ to S — 1. Hence the second-order 
energy denominators are all 2JS. 

[14] R. Shankar, Phys. Rev. Lett 45, 1088 (1980). (Our choice 
(^) corresponds to Shankar's "Hq" option.) 

[15] J. Klauder, Phys. Rev. D 19, 2349 (1979V 

[16] A BS prescription with S — > (£+1/2) in (jSJ) and a quanti- 
zation condition $ or Mt — 27r(Z + 1/2) is empirically more 
accurate; see Q. 

[17] (a). D. Loss, D. P. DiVincenzo, and G. Grinstein, 
Phys. Rev. Lett. 69, 3232 (1992); (b). J. von Delft and 
C. L. Henley, Phys. Rev. Lett. 69, 3236 (1992). 

[18] This realizes a previously hypothetical threefold symme- 
try (compare Fig. 1 of Ref. @(b).) 



TABLE I. Eigenenergies AE for small systems 



Orbit 


Exact p 


4-spin 


Theory 




N = 16 


S = 2 


5 = 8 


To (2) 


-0.151 


-0.073 


-0.0029 


Co (1,2) 


-0.272(0.006) 


-0.137 (0.0029) 


-0.0323 (0.00065) 




N = 28 


S = 7/2 


5 = 14 


To (2) 


-3.754 


-1.919 


-1.121 


Cx (1,2) 


-4.021 (0.017) 


-2.336 (0.037) 


-1.478 (0.34) 


Co (2,1) 


-4.769 (0.0003) 


-3.158 (0.0018) 


-2.033 (0.00062) 
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FIG. 1. Unit sphere representing a spin with classical or- 
bits following contours of the Hamiltonian ([h^). Orbits Coa 
(a = x,y,z) and To± sit on the special symmetry points C a 
and Tt (indicated by dots); orbits Ci a are shown by heavy 
lines, and separatrices connect the saddle points S a ,Sb,S c . 
(The unmarked saddle points are equivalent to these by the 
redundancy mentioned in text.) (b). Loop for estimating 
tunneling amplitude tcv, the shaded spherical area is indi- 
cated. Heavy dashed portions are classically forbidden and 
contribute to Im(3>c). (c). Same but for tn type orbits. 
(This area has reversed sign since this loop runs clockwise.) 
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